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We investigate interaction effects and magnetization in triangular graphene quantum dots using 
a combination of tight-binding, Hartree-Fock and configuration interaction methods. We show that 
electronic correlations play a crucial role in determining the nature of the ground state as a function 
of filling fraction of the degenerate shell at the Fermi level. A half- rilled charge neutral shell leads 
to full spin polarization but this magnetic moment can be completely destroyed by adding a single 
electron. 



Following the progress in the fabrication of graphene [1, 
0, H, Q, [1] based devices, lower dimensional structures 
such as graphene ribbons [( i III, M and more recently 
graphene quantum dots[10|, llll, LL2|, [l3|, LL4| are attract- 
ing increasing attention due to their non-trivial electronic 
and magnetic properties. In particular, it was shown that 
when an electron is confined to a triangular atomic thick 
layer of graphene with zig-zag edges, its energy spectrum 
collapses to a shell of deg enerate states at the Fermi level 
(Dirac point) (H, [H, fl7T[l8l fl^ isolated from remaining 
states by a gap. The degeneracy is proportional to the 
edge size and can be made macroscopic. This opens up 
the possibility to design a strongly correlated electronic 
system as a function of fractional filling of the shell, in 
analogy to the fractional quantum Hall effect [20], but 
without the need for a magnetic field. 

In this work, we present new results demonstrating 
the important role of electronic correlations, beyond the 
Hubbard model [H, [fl 01 and mean- field density func- 
tional theory (DFT) [17|, LL8[ . The interactions are treated 
by a combination of DFT, tight-binding, Hartree-Fock 
and configuration interaction methods (TB-HF-CI). We 
show that a half-filled charge neutral shell leads to full 
spin polarization of the island but this magnetic mo- 
ment is completely destroyed by the addition of a sin- 
gle electron, in analogy to the effect of skyrmions on 
the quantum Hall ferromagnet [2l|, [22|, [23|, l2J] and spin 
depolarization in electrostatically defined semiconductor 
quantum dotsjH, 0, 0, [HI- The depolarization of the 
ground state is predicted to result in blocking of current 
through a g raphene quantum dot due to spin blockade 
(SB)[23,[28[. 

In order to describe the interaction effects of electrons 
occupying the edge states (see Fig.l), we first solve the 
mean-field problem using a combination of tight-binding 
approach with a self-consistent Hartree-Fock method for 
the system with empty edge states, i.e. with N re f = 
N s ite — N e d ge electrons, where N site is the number of 
atoms, and N e d ge is the number of edge states. After 
extensive algebra, we arrive at the HF Hamiltonian given 
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where the operator c\ a creates a p z electron on site "i" 
with spin a. The tight-binding parameters ru a are fixed 
to their bulk values -2.5eV for nearest neighbours, and 
-O.leV for next nearest neighbours. Within the Hartree- 
Fock approximation, electrons interact with each other 
via the density matrix p = p^ +p^ from which we subtract 
the graphene bulk density matrix p bulk already present in 
the tight-binding term ru a . In addition to the on-site in- 
teraction term, all scattering and exchange terms within 
next nearest neighbours, and all direct interaction terms 
are included in the two-body Coulomb matrix elements 
(ij\V\kl) computed using Slater p z orbitals. In order to 
first empty the degenerate shell and later charge it with 
successive electrons, we transfer electrons to a metallic 
gate. The electrons in the graphene island interact with 
the gate via the term v g u given by 
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where {x^yi) are the coordinates of the atoms. This 
model assumes that the induced charge qi n d = —N e d ge is 
smeared out at positions (xi,yi) at a distance dg ate from 
the quantum dot. 

After diagonalizing the Hamiltonian,Eq.(l), we obtain 
TB+HF quasi-particles denoted by the creation operator 
b^ p , with eigenvalues e p and eigenfunctions \p > shown in 
Fig.l. We then start filling the edge states one by one, 
and solve the configuration interaction problem for the 
N a dd added electrons using 
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FIG. 1: (Color online) (a) Electronic density in a triangular 
graphene island of 97 carbon atoms where 7 electrons were 
moved to the metallic gate at a distance of d ga te. (b) Sin- 
gle particle spectrum of the structure in (a), obtained by 
tight-binding (TB, blue lines) and self-consistent Hartree- 
Fock (TB+HF, black lines) methods. The 7 edge states near 
the Fermi level are compared to DFT results. In Hartree- 
Fock and DFT calculations 7 electrons were removed, leaving 
the zero-energy states empty. The dielectric constant k is 
set to 6. Inset compares the structure of corner and side 
states obtained using Hartree-Fock and DFT calculations. In 
DFT calculations, hydrogen atoms were attached to dangling 
bonds. 
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where the indices (p, q, r, s) run over edge states, while 
the index p f runs over valence states (below the edge 
states). 

Fig. la shows the electronic density of a zig-zag edged 
triangular island of TV = 97 carbon atoms, separated 
by a distance dg ate from a metallic gate. At zero ap- 
plied voltage the island is charge neutral while applied 
voltage leads to removal/addition of electrons to the is- 
land. The single-particle energy spectrum obtained us- 
ing nearest neighbour tight-binding method (TB, blue 



lines) and a combination of the next-nearest neighbour 
tight-binding and Hartree-Fock methods (TB+HF, black 
lines) is shown in Fig. lb. As was previously shown 
by Ezawa[16], Fernandez-Rossier and Palacios[17] and 
Wang, Meng and Kaxiras|l8j using nearest-neighbour TB 
method and ab-initio DFT calculations, the linear spec- 
trum of Dirac electrons in bulk graphene collapses to a 
shell of degenerate levels at the Fermi energy, well sepa- 
rated in energy from the valence and conduction bands. 
Similar to edge states in graphene ribbons 0, 0, H, Hj], 
zigzag edge breaks the symmetry between the two sub- 
lattices of the honeycomb lattice, behaving like a defect. 
Therefore, electronic states localized on the zigzag edges 
appear with energy in the vicinity of the Fermi level. For 
a TV = 97 atoms island (see Fig.l) there are N e d ge = 7 
edge states. For the charge neutral system there is one 
electron per each edge state. A non-trivial question ad- 
dressed here is the specific spin and orbital configura- 
tion of the electrons as a function of the size and the 
fractional filling of the degenerate shell of edge states. 
Due to the strong degeneracy, many-body effects can be 
expected to be important as in the fractional quantum 
Hall effect. Previous calculations based on the Hub- 
bard approximation [TBI IhSL and local spin density 
functional theory [13, [3 showed that the neutral system 
(half-filling) has its edge states polarized. 

In order to study many-body effects within the charged 
shell of the edge states via the configuration interaction 
method, we first perform a Hartree-Fock calculation for 
the charged system of N — N e d ge electrons, with empty 
edge states and N e d ge electrons transferred to the gate, 
shown in Fig. la. The spectrum of HF quasi-particles 
is shown with black lines in Fig. lb. Due to the mean- 
field interaction with the valence electrons and charged 
gate, the degeneracy of the edge states is lifted, and the 
electron-hole symmetry is broken. The main effect is 
that a group of three states is now separated from the 
rest by a gap of ~ 0.2 eV. The three states correspond 
to HF quasiparticles localised in the three corners of the 
triangle. The same physics occurs in density functional 
calculation within local density approximation (LDA), 
shown with red lines in Fig. lb. Hence we see that the 
shell of edge states with a well denned gap separating 
them from the valence and conduction bands exists in 
the three approaches. 

The wave functions corresponding to the band of 
nearly-degenerate edge states obtained from TB-HF cal- 
culations are used as a basis set in our configuration inter- 
action calculations where we add N a dd electrons from the 
gate to the shell of degenerate states. In Fig. 2, total spin 
S of the ground state as a function of the filling number 
of the edge states is shown for different sizes of quantum 
dots. Two aspects of these results are particularly inter- 
esting. For the charge neutral case (N a dd — N e d ge = 0),for 
all the island sizes studied (N e d ge = 3 — 7), the half- filled 
shell is maximally spin polarized as indicated by red ar- 
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FIG. 2: (Color online) Spin phase diagram from TB-HF-CI 
method as a function of the size of the triangular dot charac- 
terized by the number of edge states N e d ge , and the filling of 
the zero-energy states Nadd-Nedge. Charge neutral case cor- 
responds to N a dd — N e dge = 0, for which the total spin of the 
zero-energy electrons are always maximized (S = N e d ge /2, 
indicated by red arrows). On the other hand, if the quantum 
dot is charged by 1 electron (N a dd — N e d ge = 1) then the total 
spin has minimum value, i.e. S = if N a dd is even, S = 1/2 
if N a dd is odd (indicated by blue arrows). 



rows. This is in agreement with our DFT calculations 
and reproduces previous work 17, T^. However, the spin 
polarisation is extremely fragile. If we add one extra 
electron (N a dd — N e d ge = 1), magnetisation of the island 
collapses to the minimum possible value, as indicated by 
blue arrows in Fig. 2. We note that full or partial depolar- 
ization occurs for other filling numbers but we focus here 
on the spin depolarisation at half-filling. The spin de- 
polarisation was found to be insensitive to the screening 
of electron-electron interactions, as values of dielectric 
screening constant k between 2 and 8 led to the same be- 
haviour of the charge neutral and single electron charged 
cases (not shown). 

In order to illuminate the depolarization process as an 
electron is added to the charge neutral maximally spin 
polarized system, in Fig. 3a, b we show the orbital occu- 
pancy of up-spin edge states at N a dd — N e d ge — 1, for 
the fully polarized state 5 = 6/2 (upper panel) and for 
the ground state, 5 = 0, (lower panel) for the N s u e = 97 
atoms quantum dot with 7 edge states shown in Fig.l. 
For the large spin 5 = 6/2 case, the added spin up elec- 
tron simply occupies the orbital 1 and its spin is oppo- 
site to the spins of other 7 electrons. However, the true 
ground state has 5 = 0, with spin occupancy shown in 
the lower panel. The added electron causes electrons al- 
ready present to partially flip their spin, with spin up 
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FIG. 3: (Color online) Left panel: Orbital occupancy of the 
7 edge states by spin up electrons, for the charged (N a dd — 
N e dge = 1) system, for (a) S = 6/2 and (b) S = total 
spin states. The ground state is S = (see Fig. 2). Right 
panel: Corresponding spin up-up pair-correlation functions 
(pi(ro)/>|(r)). The fixed spin- up electron is represented by a 
red arrow, and its position rO by a red circle. 



density being delocalized over all the 7 orbitals in anal- 
ogy to skyrmion-like excitations in quantum dots and 
quantum Hall ferromagnets[2l[ [22|, O, EH- The corre- 
lated nature of the 5 = spin depolarised ground state 
is illustrated by the up-up spin pair correlation functions 
given by (pf (ro)pj(r)), shown in Fig. 3c, d. The location 
of the fixed up-spin electron at site i*o is schematically 
shown with an up arrow. The up-up spin correlation 
function for 5 = 6/2 spin polarized system is strictly 
zero as there are no other spin up electrons. The spin 
correlation function for the spin depolarized ground state 
with 5 = shows exchange hole at ro , which extends to 
the nearest neighbours, and, more interestingly, for larger 
| i*o — r| , spin pair correlation function reveal a spin tex- 
ture: Beyond the exchange hole there is the formation 
of an electronic cloud with positive magnetization which 
decreases and changes sign at even larger distance, again 
consistent with the skyrmion picture [2l|, |22[. 

Experimentally, spin properties of quantum dots 
can be probed using Coulomb and spin blockade 
spectroscopy[28j. By connecting graphene quantum dot 
to leads and measuring the conductance as a function 
of gate voltage, one obtains a series of Coulomb block- 
ade peaks. The relative position of these peaks and their 
height reveal information about the electronic properties 
of the system as the number of electrons is increased. 
The amplitude of the Coulomb blockade peak is given by 
the conductivity G{ of the graphene quantum dot con- 
nected to leads via atom "i" [29j as shown schematically 
in Fig. 4a. Spin and correlation effects are reflected in 
the weight of Coulomb blockade peak proportional to the 
matrix element \(N + 1, J', S f \c\ a \N, J, 5)| 2 which gives 
the transition probability from state (N, J, 5) to state 
(TV + 1, J', 5') when additional electron is added to the 
site "i" of the graphene quantum dot from the lead. The 
ground state configuration (TV, J, 5) is controlled by the 
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FIG. 4: (Color online) (a) Schematic representation of the 
graphene island connected to the leads through a side site, 
(b) Conductivity as a function of applied gate voltage V g due 
to an electron entering the dot through the side atom, (c) 
Same as (b) but without the site dependence of the incoming 
electron. The oscillations of the spectral weight in (c) are 
purely to due correlation effects and spin blockade. 



gate voltage. For our model graphene quantum dot with 
N = 7 degenerate zero energy states, we can add a total 
of 14 electrons. Hence, one expects to obtain 14 peaks. 
In Fig. 4b some of the peaks have zero height due to spin 
blockade phenomenon. For instance, the transitions from 
(N = 7, S = 7/2) states to (N = 8, S = 0) states are 
spin blocked since it is not possible to change the spin of 
the system by AS = —7/2 by adding one electron with 
S = 1/2. Similarly, transitions from (N — 9, S — 1/2) 
states to (N = 10, S = 4/2) states are spin blocked. 
Besides the spin blockade, one sees a strongly oscillating 
feature of the spectral function heights. This is due to (a) 
strongly correlated nature of the states \N,S) , and (b) 
specific choice of the site "i" , where the lead is attached 
to. Here, we chose a site close to the middle of one of the 



sides of the triangle. The overlap of the site wavefunction 
is strongly dependent on the nature of the edge-states. In 
particular, the existence of corner states, as discussed in 
Fig.l, strongly affects the transition probabilities. To 
isolate the effect of correlation to the lead's position, in 
Fig. 4c we plot the conductivity assuming that the weight 
of the site "i" is a constant independent of the edge-state. 
As a result, the weights of spectral peaks are different, 
except for N = 8, 9, 10 where the spin blockade occurs. 
These results show how to detect the spin depolarization 
in transport experiments. Ultimately, we show here that 
one can design a strongly correlated electron system in 
carbon based material whose magnetic properties can be 
controlled by applied gate voltage. 
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